NPS  ARCHIVE 
1962 

GLUSE,  M. 


"•tavsS? 


is&ES 


wjsi 


ACCELERATION  OF  AN  UNBALANCED  ROTOR 


THROUGH  ITS  CRITICAL  SPEEDS 


1c  *k 


Michael  R0  Gluse 


ACCP-LE  RATION  OF  AS  UNBALANCED  ROTOR 


THROUGH  ITS  CRITICAL  SPEEDS 


by 

Michael  R.  Gluse 
/< 

Lieutenant ,  United  States  N&vy 


Submitted  in  partial  fulfillment  of 
the  requirements  for  the  degree  of 

MASTER  OF  SCIENCE 
IN 

MECHANICAL  ENGINEERING 


United  States  Naval  Postgraduate  School 
Monterey j,  California 


19  6  2 


.  OL 


ACCELERATION  .  F  AN  UNBALANCE1)  ROTOR 
THROUGH  ITS  CRITICAL  SPEEDS 
by 

Michael  R,  Gluse 


This  work  is  accepted  as  fulfilling 
the  thesis  requirements  for  the  degree  of 
MASTER  OF  SCIENCE 
IN 

MECHANICAL  ENGINEERING 
from  the 

United  States  Naval  Postgraduate  School 


ABSTRACT 


The  motion  of  an  unbalanced  rotor  during  acceleration  through  its 
critical  speed  is  studied  by  numerical  solutions  obtained  with  a  digi¬ 
tal  computer.  The  rotor  is  laterally  restrained  in  two  orthogonal 
directions  by  linear  springs  and  is  accelerated  by  a  constant  applied 
torque.  It  is  found  that,  for  a  fixed  combination  of  rotor  unbalance 
and  lateral  stiffnesses,  the  applied  torque  must  exceed  a  limiting 
value  in  order  to  accelerate  successfully  through  the  critical  speed 
region.  With  smaller  applied  torques,  after  initial  acceleration  from 
rest  the  speed  oscillates  continually  about  the  critical  speed  and  the 
lateral  excursions  grow  steadily.  The  conditions  necessary  for  a 
successful  acceleration  are  established  and  the  maximum  lateral  ex¬ 
cursions  during  successful  accelerations  are  determined.  The  effects 
of  small  amounts  of  viscous  damping  in  the  lateral  directions  are  also 
obtained.  Finally,  an  example  problem  is  solved  to  illustrate  the 
extension  of  the  solution  to  those  problems  where  applied  torque  is 
a  function  of  speed. 
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NOTATION 


c,  ~  Damping  in  x  direction 

c,c  -  Critical  damping  in  x  direction 

c2  =  Damping  in  y  direction 

czc  “  Critical  damping  in  y  direction 

e  =  Eccentricity  of  mass  center  from  center  of  rotation 

I  -  Moment  of  inertia  about  a  longitudinal  axis  through  the  center 

of  rotation 

I0  =  Moment  of  inertia  about  a  longitudinal  axis  through  the  mass 
center 

k,  =  Combined  shaft  and  bearing  stiffness  in  the  x  direction 

k2  =  Combined  shaft  and  bearing  stiffness  in  the  y  direction 

lbf  =  Pounds  force 

lbm  rs  Pounds  mass 

m  =  Mass  of  the  system 

M  =  Applied  torque 

2 

P(l)  =  Eccentricity  parameter  me  /I 
P(2)  =  Torque  parameter  hfcn/Tk, 

P(3)  =  Stiffness  ratio  k2/k, 

P(4)  =  Damping  parameter  (2c,  /c  =  (2^/0^) 

P(5)  s  Modifying  parameter  for  applied  torque 

R/e  =  Dimensionless  resultant  deflection 

$  -  Dimensionless  deflection  in  the  x  direction 

7^  =  Dimensionless  deflection  in  the  y  direction 

0  =  Angle  measured  counter  clockwise  from  x  axis  to  radial  axis 

originating  at  the  center  of  rotation  and  passing  through  the 
mass  center 

§  =  Dimensionless  time 

x,y  =  Coordinates  of  center  of  rotation 

y  =  Coordinates  of  mass  center 
o  )  -'o 
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1.  Introduction, 


The  maximum  amplitude  of  vibration  of  an  unbalanced  rotor  upon 
acceleration  through  its  critical  speed  is  of  real  importance  to  design 
engineers. 

An  analytic  solution  was  obtained  by  Lewis  [l]  for  the  linear 
single  degree  of  freedom  system  having  constant  acceleration  and  con¬ 
stant  force  amplitude,  with  and  without  damping,  Baker  [2]  obtained  a 
mathematical  machine  solution  of  a  linear  two  degree  of  freedom  system 
with  constant  acceleration,  also  with  and  without  damping,  Meuser  and 
Weibel  [3]  obtained  a  solution  on  the  mechanical  analyzer  for  a  single 
degree  of  freedom  system  having  constant  acceleration  and  linear  plus 
cubic  elasticity,  with  and  without  damping.  An  analog  computer  solution 
to  the  damped  and  undamped  linear  two  degree  of  freedom  system  was  ob¬ 
tained  by  McCann  and  Bennett  [4]  .  Here  again  a  constant  acceleration 
was  assumed,  Dornig  [ 5 J  solved  analytically  the  undamped  single  degree 
of  freedom  system  with  constant  acceleration. 

Both  Biezeno  and  Grammel  [6]  and  Baker  [2]  mention  the  possibility 
of  an  unsuccessful  acceleration;  that  is,  an  acceleration  to  the  vicinity 
of  but  not  through  the  critical  speed.  In  this  case,  the  energy  supplied 
is  absorbed  by  the  vibrations  and  damping,  if  present,  rather  than  be¬ 
ing  absorbed  by  the  rotor  itself  in  the  form  of  kinetic  energy  of  ro¬ 
tation. 

In  the  material  that  follows,  a  system  having  two  lateral  degrees 
of  freedom  is  investigated.  The  rotor  is  accelerated  by  a  constant 
applied  torque  rather  than  having  a  constant  angular  acceleration.  For 
this  system,  the  maximum  amplitudes  of  vibration  will  be  obtained  for  an 
appropriate  range  of  dimensionless  parameters.  Included  will  be  the 
damped  and  undamped  cases,  and  those  of  equal  and  unequal  spring  con¬ 
stants  in  the  x  and  y  directions.  In  addition,  the  areas  of  successful 

*  Square  brackets  refer  to  bibliography 
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and  unsuccessful  acceleration  will  be  defined*  Finallys  an  example 
will  illustrate  the  extension  of  the  solution  to  those  problems  where 
torque  is  a  function  of  speed* 
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2.  Equations  of  Motion. 


vV 

Figure  1  shows  the  idealized  system.  The  rotating  mass  m  has 
moment  of  inertia  I  about  the  longitudinal  axis  through  the  mass 
center  which  is  displaced  a  distance  e  from  the  center  of  the  support¬ 
ing  shaft  S,  Lateral  displacement  is  resisted  by  the  stiffnesses  of 
the  shaft  and  its  bearings,  and  and  the  associated  viscous  damp¬ 

ing  c^  and  c ^  A  driving  torque  M  accelerates  the  rotor. 

Figure  2  shows  the  free  body  diagram  of  the  system  and  includes 
all  real  and  D'Alembert  forces.  The  origin  of  the  coordinate  system 
is  at  point  B,  which  is  the  location  of  the  shaft  center  S  when  the 
rotor  is  undeflected. 

Summing  forces  and  moments: 


7,Fx  : 

A24- 

m  — 
dt2 

i- 

_  ftd2&  \  • 

=  IMS  + 

© 

(1) 

2F_y: 

m--J7 

JUtz 

+  Cl  it 

+  z  cut 

4-  ^?2  A 

/  Je  \»  .  „  ,  i-il’s  i 

(2) 

c!2©  ->,<d2A\  /AXH-\  /d2/*\  • 

2MS:  I0J^  =  i)  +  (3) 


Taking  the  moment  of  inertia  I  about  the  shaft  center  S  rather 
than  the  mass  center  reduces  Eq.  3  to: 


M  -  me 


(4) 


where  X  =  16  4- 

It  is  convenient  at  this  point  to  reduce  Eqs,  1  and  2  to  dimension¬ 
less  form  by  dividing  by  ek^,  and  Eq,  4  to  dimensionless  form  by  divid¬ 
ing  by  e^k^o 


Then  Eq.  1  becomes: 
^Figures  begin  on  page  23, 
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(5) 


rn  d? 
M  clt 


fd)  ,  ^  A  ;*] 


*1 

T 


Yr\ 

a 


J9)2  /dze 

nlMB+{7r> 


n 


Eq0  2  becomes 


X  w®  ‘  g  id) •  i(i)‘X [(§'-•  »•] 


and  Eq,  4  becomes: 


d29 


ez*W 


M  m 

e ’’^Pi  e  ^?i 


Let:  =  dimensionless  deflection  in  the  x  direction  =  «—■ 

=  dimensionless  deflection  in  the  y  direction  = 


dimensionless  time  =  \f -  £ 


4*J  __  ££/ilf  f  _fc_ 

jj-i  —  J  fT2  cbt  J  ~  ^  m  1 


dtz  ~  w 

tion  with  respect  to  dimensionless  time  §  . 
Then  Eq0  5  becomes: 


where  the  dot  denotes 


r  + 


c, 


a[  fh'ki 


y*  -h  yT  —  ( Q )  Q  -A  ©  fiMA*  Q 


Eq,  6  becomes 


ij  -h  y==r  i  ^  ^  ^  “  fe)  6-6  Cm  9 


and  Eq*  7  becomes,  after  multiplication  by 


me 


-  Mm  )  fr\e: 

©  —  ^  1  + 


lh 


f  9  -  <9  J 


(6) 


(7) 


dif ferentia- 


(8) 


(9) 


(10) 
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Examining  th  iar  \  ixhl  '•  sr.  l  .  *  and  9,  let 

C,c  ci it  leal  damping  in  the  x  direction 
Cic  s?  critical  damping  in  the  y  direction 
Substituting,  Eq*  8  becomes: 

S  i~  ^2.  ~£^'J  S  -h  f  ~  (©)  AM/  0  -f-  Q  & 

and  since: 


ZfP*7 


=  ‘Z  4  to 


(11) 


Cxc  ^2,4Pr)^l  -  zfrP^P 

1  -  <z2J=-  a  [HEP 
2  V  Mz 

Eq.  9  becomes: 


■>j  +•  (  2.  —  (B)X^  0~0  d&X.  0 


(12) 


Taking  the  damping  ratios  in  the  x  and  y  directions  to  be  equal 
(i.e.,  c]/cic  ~  c2''c  2  ^ s  ^  *“S  aPParent  from  an  examination  of  Eqs,  10, 
11  and  12  that  four  independent  dimensionless  parameters  define  the 
problem,  namely: 


P(l)  - 

r n  e 

I 

P(2)  *  - 

Mm 

P(3)  - 

Mi. 

T, 

P(4)  »  ( 

zf 

Cic 

)~(*t 

where  is  equal  to  or 

less  than 

kr 

Performing 

Eq.  10  becomes: 

&  =  P(  2)  +  P(/)[ 

S  9 

4  4 

uxy.  &J 

(13) 


5 


Eq*  1 1  )Cco  cr 

f  4-  P(4)  *S  4-  S  ~  (©)  ©  4-  <9  ©  (14) 

and  Eq.  12  becomes: 

^  H-  P(4)  P(3)  '  -(-  p(3)  7^  0  -e  un,  9  (15) 

Equations  13,  14  and  15  are  the  dimensionless  equations  of  motion 
which  will  now  be  arranged  in  suitable  form  for  solution  by  numerical 
methods . 
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3,  Numerical  Solution. 


The  Runge-Kutta  fourth  order  method  of  numerical  integration  was 

selected  to  solve  the  equations  of  motion.  This  choice  was  made  on  the 

basis  of  several  considerations,  namely:  it  is  applicable  to  non-linear 

differential  equations;  it  is  M sel f- staring” ,  i.e.,  only  the  functional 

values  at  a  single  previous  point  are  required  to  obtain  the  functional 

values  ahead;  it  is  relatively  simple  to  program;  finally,  since  it  is 

A 

essentially  the  Taylor  series  solution  through  terms  of  order  h  ,  it 
offers  an  acceptable  degree  of  accuracy  for  the  solution  of  the  problem. 
Before  proceding,  however,  two  additional  parameters  will  be  defined  in 
order  to  provide  the  equations  of  motion  with  greater  flexibility  and  to 
reduce  the  computation  time. 

The  torque  parameter  P(2)  is  a  constant  for  most  of  the  presenta¬ 
tion.  However,  to  provide  for  the  case  where  it  is  a  function  of  speed, 
let  P(2)  be  replaced  by  P(2)P(5),  where  P(5)  is  either  unity  for  the 
constant  torque  cases,  or  a  function  of  speed  for  the  variable  torque 
case. 

Also,  let: 


In  order  to  employ  the  Runge-Kutta  method  of  numerical  integration, 
it  is  necessary  to  reduce  the  three  equations  of  motion  to  a  system  of 
six  first  order  differential  equations  of  the  form: 


Let : 


Q 


(16) 


(17) 


(18) 
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Substituting  these  variables  and  the  parameters  P(2)P(5)  and  P(6) 
in  the  equations  of  motion,  Eq,  13  becomes: 


c<  - 


i'-i)  ?(z)  -  PC)!  J1  fiuK.  © 


L 


F 


& 


(19) 


Eql  14  becomes : 

-j~  r  (4-)  $  ^  &  4-  ck  ©  (20) 

and  Eq.  15  becomes: 

£  +  P(4)  P(fc)  4  P(b>o^  -y2^©  -oV  ©  (21) 

Equations  16  through  21  constitute  the  system  of  first  order  differ¬ 


ential  equations  which  must  be  reduced  to  the  form: 

dii  =£(>,  ,  ■']*(?)) 

Appendix  I  contains  this  rearrangement.  The  resulting  equations  are: 


9  -  oK 

i  =/> 

3  -  ^ 


(16) 

(17) 

(18) 


U  —  |  PM  PCs)  4-  P(i)P(4 )  P(U)  p>c*, 3  ©  P(l)  PCs)  yccxQ-  - 


PC 0  p(4)  ^  We  -  p(0  J  0 


(22) 


PC \)^z0Ma9  .C*L*e  f  P(l)  p(4)  P  /oM  9  4-  (23) 


PCO  _T  fiM-  0  Ukl  Q  P  U*~  om\.  Q  -  P( 4.)  p[u)  —  P(  3)  — 

Pto  p(s)c^l©  -  P0)^z  a  +  p(i)p(4)p(6)  p  /oi Mza  + 

PCOp^tj  ^-Qjj  £  i-pcoj 


8 


(24) 


41'  py<^  Q  -  P(4)  j'1  —  S  —  PC  l)  al1  A.U'L'Q  Caa  &  + 


P(l)  P(4)  P((o)  aU- L  6  Cao  Q  -f-  P(i )  P(3)  Ij  A-ua^B  Q3i  0  — 
PO)<AX  Cap3 9  f  P(i)  P(4)  ft1  Q><3 &■  +-  P(i)  fdezL2- 9  -h 


P(2)  P(£T_)  /Cua. 


p(0] 


Having  selected  the  integration  method  and  reduced  the  equations 
of  motion  to  the  proper  form  for  solution  by  this  method,  the  program 
itself  will  now  be  considered. 
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4.  Computer  Progran  Design  and  Use, 

Several  factors  governed  the  form  of  the  computer  program.  In  the 
material  that  follows,  the  more  important  of  these  are  discussed  in  de¬ 
tail,  while  those  of  lesser  importance  are  mentioned  briefly.  Details 
can  be  found  in  the  block  diagram  or  in  the  program  itself,  Appendices 
III  and  IV,  respectively. 

Failure  to  accelerate  through  the  critical  speed  region  can  re¬ 
sult  from  either  of  two  closely  related  conditions,  insufficient  driving 
torque  or  excessive  unbalance  of  the  rotor.  An  insufficient  driving  torque 
results  in  a  condition  wherein  the  angular  velocity  of  the  rotor  remains 
in  the  vicinity  of  that  corresponding  to  its  natural  frequency  of  vibration 
for  a  long  period  of  time.  In  this  case  the  inertia  forces,  though  initially 
small,  eventually  produce  large  lateral  deflections  which  in  turn  develop 
counter  torques  opposing  the  driving  torque.  An  excessive  unbalance  pro¬ 
duces  large  inertia  forces  as  a  result  of  the  greater  eccentricity  of  the 
mass  center.  The  resulting  lateral  deflections  also  produce  counter 
torques  opposing  the  applied  torque.  It  follows,  then,  that  the  success 
or  failure  of  an  attempted  acceleration  through  the  critical  speed  region 
is  dependent  upon  the  applied  torque  and  the  eccentricity  of  the  mass 
center,  other  physical  quantities  remaining  constant.  Hence  a  plot  of  the 
eccentricity  of  the  mass  center  versus  the  applied  torque  for  various  accelera¬ 
tion  attempts  will  indicate  points  of  either  successful  or  unsuccessful 
accelerations  through  the  critical  speed  regions.  A  boundary  line  may 
then  be  drawn  to  separate  those  points  representing  successful  accelerations 
from  those  points  representing  unsuccessful  accelerations. 

For  a  given  system,  the  parameter  P(l)  is  directly  proportional  to  the 
square  of  the  eccentricity  of  the  mass  center,  while  parameter  P(2)  is 
directly  proportional  to  the  applied  torque.  These  parameters,  the  eccentri¬ 
city  parameter  P(l)  and  the  torque  parameter  P(2),  are  the  ordinate  and 
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abscissa,  respectively,  of  the  curve  which  is  the  boundary  between  the 
areas  of  successful  and  unsuccessful  accelerations.  The  boundary  it¬ 


self  is  defined  by  pairs  of  points.  Each  pair  of  points  represents  two 
acceleration  attempts,  one  successful,  the  other  unsuccessful,  differ¬ 
ing  by  five  per  cent  of  the  torque  parameter  P(2),  all  other  parameters 
being  constant.  It  is  now  appropriate  to  define  a  successful  and  an  un¬ 
successful  acceleration. 

Referring  to  Eqs.  20  and  21,  one  of  the  natural  circular  frequencies 
of  vibration  of  the  dimensionless  system  is  unity  and  the  other  is  the 


The  maximum  value  of  the 


latter  is  unity.  On  this  basis,  a  successful  acceleration  was  defined 


as  an  acceleration  in  which  a  (dimensionless)  angular  velocity  of  two 


(twice  that  corresponding  to  the  higher  natural  circular  frequency  of 
vibration)  was  reached.  That  is,  if  a  dimensionless  angular  velocity  of 
two  was  reached,  the  acceleration  through  the  critical  speed  region  was 
considered  successful  and  the  program  terminated. 

The  possibility  existed,  however,  that  during  the  fluctuations  in 
the  angular  velocity  of  an  otherwise  unsuccessful  acceleration  a  speed  of 
two  might  have  been  momentarily  attained,  upon  which  the  program  would 
have  terminated  and  the  acceleration  would  have  been  erroneously  termed 
successful.  To  investigate  this  possibility,  several  random  successful 
and  unsuccessful  accelerations  were  initiated,  and  graphs  of  dimensionless 
angular  velocity  versus  dimensionless  time  were  made.  Two  such  curves  can 
be  seen  in  Fig.  3.  For  the  successful  acceleration  the  torque  parameter  is 
five  per  cent  larger  than  for  the  unsuccessful  acceleration,  all  other  para¬ 
meters  being  unchanged.  It  is  apparent  that  the  maximum  angular  velocity 
for  the  unsuccessful  acceleration  is  well  below  that  value  (2.0)  which  would 
terminate  the  program  and  falsely  indicate  a  successful  acceleration. 
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An  unsuccessful  acceleration  was  defined  essentially  in  terms  of  a 
successful  acceleration.  Starting  from  a  known  successful  acceleration 
the  torque  parameter  P(2)  was  decreased  in  five  per  cent  increments,  while 
the  time  to  reach  a  dimensionless  angular  velocity  of  two  was  recorded. 

In  regions  far  removed  from  the  boundary  between  the  successful  and  un¬ 
successful  acceleration  regions,  the  time  required  for  succeeding  success¬ 
ful  accelerations  increased  by  five  per  cent.  In  the  vicinity  of  the 
boundary,  however,  the  additional  time  required  for  a  successful  accelera¬ 
tion  exceeded  the  expected  five  per  cent,  but  in  no  case  did  it  exceed 
thirty  per  cent.  On  this  basis,  the  time  allowed  for  each  successive 
acceleration  was  twice  that  required  by  the  previous  successful  acceleration. 
If  a  dimensionless  angular  velocity  of  two  was  not  achieved  within  this 
period  of  time,  the  attempt  was  considered  unsuccessful  and  the  program 
terminated.  Figure  3  also  illustrates  this  point.  The  upper  curve  re¬ 
presents  a  successful  acceleration  in  which  a  dimensionless  angular  velocity 
of  two  was  reached  in  time  t.  The  lower  curve  represents  an  unsuccessful 
acceleration  in  which  a  dimensionless  angular  velocity  of  two  was  not 
achieved  in  time  2t.  The  torque  parameters  for  these  attempts  differ  by 
five  per  cent,  all  other  parameters  being  constant.  Hence  these  two 
attempts  correspond  to  one  of  the  several  pairs  of  points  which  define 
the  boundary  between  regions  of  successful  and  unsuccessful  accelerations. 

For  the  undamped  case,  an  energy  balance  affords  a  check  on  the  valid¬ 
ity  of  the  results.  At  any  given  time  the  total  energy  supplied  to  the 
system  must  equal  the  sum  of  the  kinetic  energies  of  rotation  and  trans¬ 
lation,  and  the  potential  energy  stored  in  the  springs,  or: 


Ma  = 


4  Y-o 
dJt 


4-  -L  rvj 
Z 


(25) 
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An  equivalent  statement  in  dimensionless  form  is. 


ru)  e  - 


P(0 


(e) 


.  (i-r  (.(ij'r  4  f 2 


n 


L  P  0 

©  (  S  b  ut,  9  -  Os.  6> 


4-  Pfs)  >r  - 

I 


(26) 


For  the  undamped  runs,  this  check  was  performed  at  the  end  of  each 
acceleration;  ice.,  employing  the  terminal  values  of  the  quantities 
in  Eq.  26. 

The  choice  of  a  suitable  time  increment  was  made  on  the  basis  of 
a  comparison  of  the  results  from  several  solutions  of  the  same  problem 
employing  various  time  increments.  Details  are  given  in  Appendix  II. 

A  (dimensionless)  time  increment  of  0.10  was  selected. 

Since  the  objectives  of  this  investigation  included  a  determination 
of  the  maximum  amplitudes  of  vibration  upon  acceleration  through  the 
critical  speed  region,  this  quantity,  the  resultant  of  the  deflections  in 
the  x  and  y  directions,  was  computed  each  time  increment  and  compared  with 
a  stored  maximum  which  had  been  determined  earlier,  the  larger  being  re¬ 
tained. 

Provision  was  also  made  to  terminate  the  solution  at  any  point  and  to 
resume  the  solution  at  this  same  point  at  a  later  date.  This  innovation 
provided  a  considerable  savings  in  computer  time  since  any  time  period, 
however  small,  could  be  completely  utilized,  and  no  solutions  were  lost 
because  of  unexpected  requests  to  release  the  computer.  The  termination 
and  resumption  was  executed  by  writing  the  values  of  the  (dimensionless) 
deflections  and  velocities  in  the  x  and  y  directions,  the  angular  velocity 
and  displacement  and  elapsed  time  on  magnetic  tape  in  binary  form,  and 
reading  these  same  values  back  into  the  program  by  the  setting  of  certain 
selective  jump  switches. 

By  way  of  a  check  of  the  computer,  the  output  tape  unit  and  the  tape 
itself,  the  same  short  run  was  performed  at  the  beginning  of  each  period 
of  computer  use. 
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The  program  output  included  the  parameters  employed,  the  results  of 
the  energy  balance,  the  maximum  angular  velocity  and  the  time  of  its 
occurrence,  ind  the  final  values  of  the  deflections  and  velocities  in  the 
x  and  y  directions,  angular  velocity,  revolutions,  and  elapsed  time.  Also 
recorded  were  the  maximum  amplitudes  of  vibration  in  the  x  direction,  y 
direction  and  the  maximum  resultant  amplitude.  For  each  of  the  three 
maximum  amplitudes  of  vibration,  the  time  of  its  occurrence,  angular  velo¬ 
city  and  the  x  direction,  y  direction  and  resultant  amplitudes  of  vibra¬ 
tion  were  also  recorded. 

In  general,  the  investigation  proceeded  as  follows.  Selecting  an 
appropriate  set  of  values  for  parameters  P(l),  P(3)  and  P(4),  the  torque 
parameter  P(2)  was  varied  to  define  the  boundary  between  the  regions  of 
successful  and  unsuccessful  accelerations.  Having  defined  the  boundary,  P(2) 
was  then  increased  through  an  appropriate  range  of  values  to  obtain  the 
corresponding  maximum  amplitudes  of  vibration.  This  procedure  was  then 
repeated  for  different  values  of  P(l),  P(3)  and  P(4)0 

The  following  section  presents  the  results  of  the  investigation. 
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5.  Discussion  of  Results. 


The  results  of  the  investigation  are  presented  as  follows. 

First,  the  boundary  between  the  regions  of  successful  and  unsuccessful 
acceleration  is  discussed  for  stiffness  ratios  of  unity,  0,50,  0.75  and 
0.25.  Then  the  curves  of  the  maximum  dimensionless  amplitudes  of  vibra¬ 
tion  are  presented.  Finally,  the  effects  of  added  damping  are  considered. 

A  rotor  having  equal  spring  constants  in  the  x  and  y  directions  (P(3) 

=  1.0)  and  no  damping  (P(4)  =  0)  was  investigated  first.  Figure 4*  indicates 
the  location  of  the  boundary  between  the  regions  of  successful  and  un¬ 
successful  accelerations  as  a  function  of  the  eccentricity  parameter  P(l) 
and  the  torque  parameter  P(2).  This  curve  is  essentially  a  reference  for 
the  boundary  curves  that  follow  since  unequal  stiffness  and  added  damping 
are  more  conveniently  discussed  in  terms  of  their  effects  on  this  particu¬ 
lar  result. 

Figure  4  shows  that  at  successively  smaller  eccentricities,  smaller 
applied  torques  are  required  for  successful  accelerations,  since  the 
smaller  resultant  inertia  forces  produce  smaller  counter  torques. 

Decreasing  the  stiffness  ratio  to  0.50,  the  undamped  boundary  curve 
was  found  to  be  as  shown  in  Fig.  5.  For  any  non- unity  stiffness  ratio, 
the  system  will  have  two  critical  speed  regions;  in  this  particular  case 
one  is  in  the  vicinity  of  a  (dimensionless)  angular  velocity  of  /f  0.501 
and  the  other  is  in  the  vicinity  of  unity.  It  follows,  then,  that  an 
unsuccessful  acceleration  can  result  from  an  inability  to  accelerate 
through  either  the  lower  critical  speed  region  or  the  higher  critical 
speed  region.  Figure  5  illustrates  both  of  these  conditions.  The  bound¬ 
ary  between  the  regions  of  successful  and  unsuccessful  accelerations  is 
seen  to  vary  between  two  parallel  limits.  In  the  vicinity  of  the  higher 
limit,  the  unsuccessful  acceleration  points  defining  the  boundary  represent 
failure  to  accelerate  through  the  higher  critical  speed  region.  Similarly, 
*Figure  4  can  be  found  on  page  26, 
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those  unsuccessful  accelerations  defining  the  boundary  in  the  vicinity  of 
the  lower  limit  represent  failure  to  accelerate  through  the  lower  critical 
speed  region. 

The  limits  mentioned  above  have  a  definite  relationship  to  each  other. 
For  a  given  eccentricity,  the  ratio  of  the  lower  limit  to  the  upper  limit 
is  very  nearly  equal  to  the  stiffness  ratio,  This  can  be  explained 

M  m 

by  examining  the  dimensionless  torque  parameter  P(2)  =  * 

I 

It  is  apparent  from  this  dimensionless  grouping  that  a  system  which  has  a 
spring  constant  of,  say,  Ak^  will  require  a  torque  of  AM  for  a  given 
eccentricity.  In  the  vicinity  of  the  lower  critical  speed  region  of  the 
system  under  consideration,  the  system  behaves  essentially  as  a  single  de¬ 
gree  of  freedom  system  having  a  stiffness  k^  =  0.50k^„  Hence  the  torque 
required  for  acceleration  is  0.50M.  Since  the  boundary  curves  are  plotted 
versus  P(2),  which  is  a  function  of  kp  the  substitution  of  0.50M  for  M 
results  in  a  decrease  of  the  torque  parameter  from  P(2)  to  0„50P(2). 

The  fact  that  a  portion  of  the  boundary  curve  is  in  the  vicinity  of 
the  lower  limit  shows  that,  in  this  region,  a  successful  transit  of  the 
lower  critical  speed  region  will  result  in  a  successful  transit  of  the 
higher  critical  speed  region  also.  This  condition  may  be  due  to  a 
favorable  phase  relationship  between  those  vibrations  induced  at  the  lower 
and  higher  critical  speeds.  That  is,  those  vibrations  remaining  in  the 
system  as  a  result  of  the  transit  of  the  lower  critical  speed  region  exert  a 
forward  torque  opposing  that  developed  by  the  vibrations  induced  in  the 
higher  critical  speed  region  to  the  extent  that  a  smaller  torque  parameter 
suffices  for  a  successful  acceleration. 

For  a  stiffness  ratio  of  0.75,  the  undamped  boundary  curve  was  found 
to  be  as  shown  in  Fig.  6.  The  general  shape  of  the  curve  is  comparable  to 
that  of  the  stiffness  ratio  previously  discussed,  0.50;  however,  some  dif¬ 
ferences  do  exist. 
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The  boundary  is  seen  to  vary  between  two  limits  much  the  same  as 
the  boundary  of  the  0.50  stiffness  ratio  case;  however,  for  a  given  eccentri¬ 
city,  the  ratio  of  the  torque  at  the  lower  limit  to  that  at  the  upper  limit 
is  influenced  considerably  by  the  close  proximity  of  the  two  natural  fre¬ 
quencies  of  vibration.  In  the  higher  torque  regions,  this  influence  is 
quite  apparent,  and  as  a  result,  the  limit  ratio  differs  somewhat  from  the 
stiffness  ratio.  In  the  lower  torque  regions,  agreement  between  the  stiff¬ 
ness  ratio  and  the  limit  ratio  is  quite  good  since  there  are  many  revolu¬ 
tions  of  the  rotor  separating  the  critical  speed  regions  and  their  inter¬ 
action  is  reduced  considerably. 

For  a  stiffness  ratio  of  0.25,  the  critical  speeds  of  0.50  and  unity 
are  sufficiently  displaced  from  each  other  that  little  or  no  interaction 
of  vibrations  is  experienced.  The  counter  torques  developed  in  the  higher 
critical  speed  region  are  larger  than  those  induced  in  the  lower  critical 
speed  region  due  to  the  larger  angular  velocity.  It  follows  that  an  un¬ 
successful  acceleration ’ represents  a  failure  to  accelerate  through  the 
higher  critical  speed  region.  The  unsuccessful  acceleration  points  defining 
the  boundary  curve  shown  in  Fig.  7  represent  just  such  a  failure. 

By  way  of  a  brief  summary,  the  undamped  boundary  curves  described 
above,  i.e.,  for  stiffness  ratios  of  1.00,  0.50,  0.75  and  0.25,  are 
plotted  on  a  single  graph,  Fig.  8. 

Except  as  noted  below,  the  starting  phase  0=0  was  used.  When  the 
lateral  stiffnesses  are  equal  (P(3)  =  1),  the  resulting  circular  symmetry 
assures  that  neither  the  location  of  the  boundary  nor  the  maximum  vibra¬ 
tion  amplitude  will  be  affected  by  starting  phase.  Even  with  unequal 
stiffness  ratios,  it  is  reasonable  to  expect  that  the  effect  of  starting 
phase  is  negligible  if  the  rotor  completes  a  large  number  of  revolutions 
before  reaching  the  lower  critical  speed.  This  condition  is  met  for 
small  values  of  the  torque  parameter. 
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An  anomalous  behavior  in  the  higher  torque  region  did  result  in  a 

limited  study  of  the  effect  of  starting  phase.  With  a  stiffness  ratio 

-3 

of  0,50,  it  was  found  at  P(l)  =  2o30xl0  that  the  minimum  value  of  P ( 2 ) 

for  successful  acceleration  was  46,0  per  cent  less  than  at  immediately 

adjacent  higher  and  lower  values  of  P(l).  Further  runs  were  made  at 
-3 

P ( 1 )  =  2.30x10  with  starting  values  of  0  =  45*  90  and  135  degrees  and 

it  was  found  that  the  minimum  P(2)  for  successful  acceleration  approached 

the  results  obtained  for  neighboring  values  of  P(l),  Because  the  principal 

engineering  interest  in  these  results  is  in  selecting  a  combination  of 

parameters  that  will  assure  a  successful  acce lerat ion*  regardless  of  start- 

-3 

ing  phase,  the  "stray"  point  at  P ( 1 )  =  2,30x10  obtained  with  9=0  was 
ignored  in  plotting  the  boundary  curve  of  Fig,  5. 

Figures  9  through  12  are  the  undamped  maximum  amplitude  curves  for 
the  stiffness  ratios  investigated.  These  curves  indicate  that,  for  a 
given  eccentricity,  successively  larger  applied  torques  result  in  de¬ 
creasing  amplitudes  of  vibration.  These  smaller  amplitudes  are  to  be  ex¬ 
pected  since  the  larger  torques  result  in  a  shorter  period  of  time  in  the 
critical  speed  regions. 

These  curves  further  show  that  in  regions  of  successful  accelerations 
removed  from  the  boundary  between  the  regions  of  successful  and  unsuccessful 
accelerations,  the  maximum  dimensionless  amplitudes  of  vibration  are  in¬ 
dependent  of  the  eccentricity  parameter  P(l) ;  or,  for  a  given  torque  para¬ 
meter  P(2)  : 

R/e  =  Constant 

from  which: 

R  =  (e)  (Constant) 

Therefore,  the  maximum  amplitude  of  vibration  is  directly  proportional  to 
the  eccentricity  e  of  the  mass  center.  This  is  to  be  expected  since  these 
vibrations  are  a  consequence  of  the  inertia  forces  which,  in  turn  are  direct¬ 
ly  proportional  to  the  eccentricity  of  the  mass  center. 
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The  unusual  variations  in  the  amplitude  curves  of  Figs  10  and  11 
are  a  result  of  the  transition  of  the  boundary  curve  from  one  of  its 
limits  to  the  other  limit* 

Figures  13  through  20  indicate  the  effects  of  the  addition  of  one 
per  cent  of  critical  damping  on  the  boundary  and  amplitude  curves*  With 
one  exception,  shown  in  Fig*  13,  these  curves  reveal  a  decrease  in  the 
maximum  amplitudes  and  a  decrease  in  the  torque  required  for  a  successful 
acceleration.  The  decreased  amplitudes  result  in  smaller  counter  torques, 
hence  acceleration  through  the  critical  speed  region  is  achieved  with  a 
smaller  applied  torque*  The  larger  shift  of  the  boundary  in  the  lower 
torque  regions  is  a  result  of  the  larger  number  of  cycles  during  which 
the  damping  has  acted* 

Figure  15,  however,  shows  that  the  effects  of  added  damping  are  un¬ 
favorable  for  a  considerable  range  of  values  of  eccentricity  and  torque 
parameters.  It  is  possible  that  the  added  damping,  though  reducing  the 
individual  amplitudes  of  vibration  induced  in  the  lower  and  higher  critical 
speed  regions,  could  have  also  resulted  in  a  particularly  unfavorable  phase 
relationship  between  these  vibrations*  As  a  consequence  the  resulting 
amplitudes  and  counter  torques  could  have  been  particularly  large  and  a 
larger  applied  torque  would  be  required  for  successful  acceleration* 

As  was  the  case  with  Figs*  10  and  11,  the  unusual  variations  of  the 
amplitude  curves  of  Figs.  18  and  19  are  a  result  of  the  transition  of  the 
boundary  curve  from  one  of  its  limits  to  the  other* 
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6,  Sample  Problem, 

The  computer  program  is  not  limited  to  accelerations  with  constant 
applied  torques0  The  following  is  a  sample  problem  illustrating  the  ex¬ 
tension  of  the  solution  to  a  system  having  a  variable*  rather  than  a 
constant,  applied  torque. 

An  unbalanced  rotor*  its  supporting  shaft  and  drive  sheave  have  a 

combined  mass  m  =  12,0  lbm,  and  a  moment  of  inertia  about  the  mass  center 
2 

I  =  48,0  lbm-in  .  The  mass  center  of  the  system  is  displaced  a  distance 
e  =  0.02  in. from  the  center  of  rotation.  The  stiffnesses  of  the  shaft  and 
its  supports  are  the  same  in  the  x  and  y  directions  and  are  =  100 

lbf/in.  The  drive  motor  for  this  system  has  a  starting  torque  M  =  40.0  lbf- 
in*  which  decreases  linearly  with  angular  velocity  to  a  value  of  zero  at 
55.2  rpm,  or  twice  the  critical  speed  of  the  system.  If  damping  is  negli¬ 
gible,  determine  the  maximum  amplitude  of  vibration. 

The  axis  of  the  moment  of  inertia  must  be  translated  from  the  mass 
center  to  the  center  of  rotation.  Since: 

T  =  To  +  ^e2"  , 

t-  ( 48.o  Ikm  -  6u)  i  (12.0  lb^)(o.02  =  47,0  ibn.j*.* 


Defining  the  parameters: 

P(0=  npel 

r 


Similarly : 


pc 0  =  (12.0  -Q2 

( 4.8.0  Ibm 


—  1-00  X  i O 


-4 


P(z)  - 


M  m 

TMi 


—  0.  too 


P(2)  -  (40.0  lbm) 

L48.0  joo 

This  value  of  torque  decreases  linearly  to  zero  at  a  (dimensionless)  angular 


velocity  of  2.00  (twice  critical  speed)*  therefore  the  modifying  parameter 
P(5),  which  has  been  unity  for  the  constant  torque  case,  now  becomes: 


(,> = ( 


p<-~ 


JL 

z 


and  the  dimensionless  torque  terms  in  the  equations  of  motion  become: 
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Also  : 


p(i)  !  (?)  =.  0.1  00(  I-  -§) 


■kz 

P(3)  =  -j-  -1.00 

and,  since  damping  is  negligible: 

P(4^  =  0.00 

Inserting  these  parameters  into  the  computer  program,  the  resulting  solution 
indicates : 


R 


=  lAci 


from  which: 


R - O,  j  (cjO  /yVU. 


where  R  is  the  maximum  resultant  amplitude  of  vibration. 
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Fig,  1  Idealized  system 
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Fig-  2  Free  body  diagram 
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Dimensionless  Angular  Velocity 
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Fig,  3 


Plot  of  dimensionless  angular  velocity  versus  dimensionless  time  for  a 
successful  and  an  unsuccessful  acceleration.  The  successful  acceleration 
has  a  torque  parameter  five,  per  cent  larger  than  that  of  the  unsuccessful 
acceleration,  all  other  parameters  being  unchanged. 


Figo  4  Location  of  the  boundary  between  the  regions  of 

successful  and  unsuccessful  accelerations,  P(3)  = 
1.0  (unity  stiffness  ratio)^  P(4)  »  0  (undamped 
case) o 
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Fig.  5  Location  of  the  boundary  between  the  regions  of 
successful  and  unsuccessful  accelerations  for  a 
stiffness  ratio  P(3)  -  0.50,  P(4)  ~  0  (undamped 
ease ) . 
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stiffness  ratio  P(3)  =  0.75;,  P(4)  =  0  (undamped 
case) o 
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stiffness  ratio  F(3)  -  Oc.25,  P{4)  -  0  (amdamped 
case). 
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Fig.  8  Location  of  the  boundary  between  the  regions  of  successful 
and  unsuccessful  accelerations  for  stiffness  ratios  of 
1.00,  0.75,  0.50  and  0.25,  P(4)  s  0  (undamped  cases). 
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Maximum  dimensionless  amplitude  of  vibration 
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Fig.  9  Maximum  dimensionless  amplitudes  of  vibration  as  functions  of 

eccentricity  parameter  P ( 1 )  and  torque  parameter  P{2);  stiffness 
ratio  P(3)  =  1.0*  P(4)  -  0  (undamped  case). 


Maximum  Dimensionless  Amplitude  of  Vibration 


Fig.  10  Maximum  dimensionless  amplitudes  of  vibration  as  functions  of  eccentricity  para¬ 
meter  P(l)  and  torque  parameter  P(2);  stiffness  ratio  P(3)  -  0S50*  P(4)  =  0  (un¬ 
damped  case). 


Maximum  Dimensionless  Amplitude  of  Vibration 


Fige  11  Maximum  dimensionless  amplitudes  of  vibration  as  functions  of  eccentricity  parameter 
P(l)  and  torque  parameter  P(2),  stiffness  ratio  P(3)  *  0,75*  P(4)  a  0  (undamped 
case). 
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Fig.  13  Effect  of  1„00  per  cent  of  critical  damping  on  the  location 
of  the  boundary  between  regions  of  successful  and  unsuccess¬ 
ful  accelerations,  stiffness  ratio  P(3)  -  1.00,  P(4)  »  0.02, 
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Fig.  14  Effect  of  1*00  per  cent  of  critical  damping  on  the 
boundary  between  the  regions  of  successful  and  un¬ 
successful  accelerat ions*  stiffness  ratio  P(3) 
0,50,  P (4)  »  0.02, 
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between  the  regions  of  successful  and  unsuccessful  accelera¬ 
tions,  stiffness  ratio  P(3)  =  0.75,  P(4)  =  0.02. 
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tions,  stiffness  fatio  ?(3)  -  0.25,  P{4)  -  0.02, 
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Maximum  Dimensionless  Amplitude  of  Vibration 


Fig.  1?  Maximum  di  lensionless  amplitudes  of  vibration  as  functions  of  eccentricity 

parameter  P(l)  and  torque  parameter  P(2),  P(3)  ~  1.0  (unity  stiffness  ratio). 
1,00  per  cent  critical  damping  added  (P(4)  -  0.02)„ 


Maximum  Dimensionless  Amplitudes  of  Vibration 


Fig,  18  Maximum  dimensionless  amplitudes  of  vibration  as  functions  of  eccentricity  parameter  P(l, 
and  torque  parameter  P(2)s  stiffness  ratio  P(3)  •=  0.50}  1.00  per  cent  critical  damping 
added  (P('t)  -  0.02), 
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Ftge  19  Maximum  dimensionless  amplitudes  of  vibration  as  functions  of  eccentricity  parameter  P^l 
and  torque  parameter  P(2)s  stiffness  ratio  P(3)  =  0.75.  1.00  per  cent  of  critical  damp 
ing  added  (P(4)  -  0.02). 


Maximum  Dimensionless  Amplitude  of  Vibration 
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CONVERSION  OF  EQUATIONS  FOR  NUMERICAL  INTEGRATION 
The  system  of  first  order  equations  to  be  rearranged  is  restated 
for  convenience. 


—  oL 


(16) 


^  p 

i  =  j1 


(17) 


(18) 


k  =■  PO-)  P(S)  +  PCI)  r  ^  Al 


B  -  p>  f&A  & 


(19) 


^  -f-  P(4)  tP  -f  J  —  CJr^  Q  -f-  cA  AU<~  Q 


(20) 


fS  4*  p(4>  p(fe)  £  +-  P(3)'r^  =  o< 1  0  ©  (21) 

Equations  16,  17  and  18  have  the  required  form: 

jA  =  f  (*>  t' ^00) 

It  remains,  then, to  reduce  Eqs.  19,  20  and  21  to  this  same  form0 
Rearranging  these  equations: 

<k  ~h  PC  0  B  ^  PO)  Q  ft  —  A 


©  cC  -f- 


=  B 


A 

r 


=  c 


where: 


A  - 

P(2)  PCs) 

B  = 

A2  /X-  0 

p(4)  ^  -  / 

C  = 

1  & 

-  PCs)  p  -  P(A)P(io)  p 
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Employing  Cramer's  Rule, 


c 


n 


I  PO)  f-frd  Q 

O-lM.  o  O 

f  /-a  ©  I 


—  f ( r  — 


0 


-  PC  I  )  fiAM,  Q 


o 


Thus  : 


o(.  ~ 


C  PO)  Aa-'Xj  0  ~  B  p(0  i g-ua  ©  —  A 
PCOctJe  i  po)^e  -i 


—  ©  -  P(3)t|  -  P(4)P(C>)  jB~j  P(t)  d&d  &  — 

jC-Z'r  0>  ©  —  PLA)  ^  -  JT~|  P(l)  AUL  Q  -  p(2)  P(5) 


or : 


^  = 


PO )- 


p(0  -J1  <4Ua  &  C&d.  &  -  PL  Pi  Pi 3)  7j  O  - 

P(I)  P(4)PCb)  p  G  -  P'(l)^z  aUu  ©  ©  -f- 

P(l)  H4)  ^  A^o  ©  4  PO)  J  />uc  ©  - 


PC*)  PCs) 


PCD  -I 


] 


Finally: 
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>  =  r  Ri)P(»  r  +  P(/)P(4)  P(b)  (5  ©  ~  (22) 

ft.  A  K4)  ,H  ,cUa  9  -  P(i)  S  A  A  PRft  P£r) 

r  i  -  p(o] 

L  J 


I  A 

—  ama-  ©  B 

ft*!-  &  C 

[pco-n 

u'  _i 


—  P(/)  6 


O 


Thus : 


3  -  |  A  cm  ©  +  C  PCI)  aua  z9  -h  5  PCO  Am,  ©  A^  ©  — 


or : 


3  -{  P(2-)  PftO  A^e  &  -  P(3)  7|  -  P(4)P(» 

P0)A^Z©  -P(4)(T  -I 


^U4  $ 


■A  ^  e  +  P(3)  >i  +  P(4)  «t)  pj  j  |>co-/J 


Finally : 


p=|j_P(2')PC S')  C&L  Q  —  PCOft  L/yuA  ~  &  -f  P(0  P(3)  ~yj  /Qma.2  Q+.  (23) 

P(0  P(4)  PU)  p  A^A28  ~  P(/)AZ  AH  0  A 

PC  0  P(4)  ^  A^  0  pw  0  -f  P(0  /  xi^m-  &  /C&a  Q  A 
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4  -  ' 


k=>)  >,  -  p(a)  KO  /fil  /[  I 
v  j  y'  L 


p(0 


_J 


A!  i,v  0 


P(l)  CM;  © 


O 


K1  = 


I  c^c,  0 


A 

B 

C 


[  PO)  -  I  J 


Thus : 

ii  =  6 P(i)  ^2e  -  A  w.  6  -  S  -  c  PO)  Aw.  0  to  ej  / 
[  PCli-ll 

or: 

AtfO  0  -  RO  i1  -  tJ  PO)  -  PC 2)  PCs-)  Q 

-p  P(4)d*  4-  f  4-^^/CW.©  ~p( s)l|  - 

P(4)  P((o)  p 

Finally : 

ii  =-  po)°^^l3£  p  PC)  P(4^  c^&  f  P0)j>^e  + 

P(2_ )  PCS )  <9  -f~  cK*~  j&fi'X  &  ^  P ~~  S*  ~~  P(t)  * 

^  Q  C-Uj  Q  4-  PO)  P^5)  ~7j  fl'tiA  Q  /C^hl  ©  4- 

P(/)  P{4)  P<%)  j?j  fiUA-  Q  /CJrx  Q 


J 


PO)  /Hil  0  /Uk2  & 


PO)  -ll 


(24) 
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APPENDIX 
SEILCr:ON  OF  TIME  INCREMENT 

The  selection  of  the  tine  increment  h  was  made  by  selecting  several 
representative  sets  of  parameters  and  obtaining  the  corresponding  solu¬ 
tions  for  various  h.  This  afforded  two  comparisons*  one  of  the  energy 
balances  and  another  of  the  resultant  deflections*  Then  the  torque  para¬ 
meter  was  varied  until  the  boundary  between  the  regions  of  successful  and 
unsuccessful  accelerations  was  defined*  This  afforded  a  third  check*  Re¬ 
presentative  results  for  one  set  of  parameters  are  presented  in  the  accompany¬ 
ing  table*  The  time  increment  of  0.1  was  selected  as  the  reference  since 
smaller  time  increments  produced  no  change  in  any  of  the  three  comparisons 
made  for  each  set  of  parameters. 


Per  cent  difference  as  compared 
with  results  for  h  -  0.1 

h 

Resultant 

Energy 

Boundary 

0,1 

0,00 

0.00 

0.00 

0.2 

0.00 

0,00 

5.00 

0,  3 

0.61 

0,00 

0,00 

0,4 

2.23 

2.40 

0.00 

0.5 

1,82 

7.20 

0.00 

Effect  of  time  increment  h  on  resultant 
deflection,  energy  balance  and  location 
of  the  boundary  between  regions  of 
successful  and  unsuccessful  accelerations* 


4? 


ppe:._  x  :r 

BLOCK 


Start 


u 


bet  parameter 
t  for  test  run 


u 


A 


OH  Set  initial  values 


i - 

I  Evaluate  integrals 


Compute  resultant 
deflection  R 
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Mo 


Ut 


Print  desired 
outputs 


Set  new 
parameters 


T~> 


5- 


PRO OR AM  B2 

DIMENSION  Y  (  6  )  ,  0(6),  P(6),  Y  F  { t> ) 

COMMON  P 

UUN=0. 

Z  =  0. 

(>{  1  )  =  .0  1 
P( 2  )=.07 
P  (  3  )  =  I  . 

P(  4  )=(;. 

R=0. 

DELP=0. 

1-1=  „  1 

RUN  =  R«JN+  I  . 

P( 5  )  =  1  . 

P(6)  =  S«RTHP(4)  ) 

SP  0 1-,  X  =  0  . 

XMSPO=0. 

Y (  1  )  =  U. 

Y( ?  )=0. 

Y  (  3  )  =  0  . 

Y  (  4  )  -  0  . 

Y( 5 )=0. 

Y  (  6  )  =  0  . 

X=0  - 

Yt (  i  )=0. 

Y  Z  (  2  )  =  0  o 

Y  ft  (  3  )  =  0  * 

CALL  RKUTT A ( 6  »  X  »  Y , H ) 

X-X  +  H 

R=  SORT  F ( Y (  1  >*Y(  1  )+Y(2) »Y(?) ) 

I  F ( SPCMX— Y ( 6 )  )  75,20,20 
SPDNX=Y( 6 ) 

XMS  PD=  X 
RFV=Y( 31/6.28 

1  F  (  A  B  S  F  (  Y  (  1  )  J-ARSF (YE(  I  )  )  )  \  ,  1  ,4 
IF(ABSF(Y(2)  )-AH5F(YZ(?)  )  12,2,5 
IF  (  A  B  S  F  (  R  )  -  A  B  S  F  {  Y  R  (  3  )  )  )  5 ,3, 6 
YL (  I  )=Y(  i  ) 

YF ( 2 )-Y( 2  ) 

Yt-  (  3  )  - R 
YF ( 4 ) =  Y ( 6 ) 

Y  t ( b )=X 
GO  TO  1 


Y  Z ( 5 )  ,  YR ( 5 ) 


COMPUTER  PRC  1  RAM 


YZ ( 2 )=Y( 2) 

Y  Z  (  3  )  =  R 

YZ  (  4  )=Yt  6) 

YZ(S)=X 
GO  TO  2 

6  YR{  l)=Yt  1  ) 

Y  R  (  2  )  -  Y  {  2  ) 

YR{ 3 ) =R 

Y« ( 4  )  -  Y  (  6  ) 

YR ( S )=X 

5  IF  (SENSE  SWITCH  2)  9  I  ,  92 

91  PAUSE  2 

REAU  TAPE  2, ( Y( 1 ) , 1= 1 , 6) ,X 
REWIND  2 

92  IE (SENSE  SWITCH  1)  22 ,7 

!  IF (Y( 61-2. ) 10,8,8 

22  PRINT  23 

23  FORE  AT l  I6HTIME  CUTOFF  USE  T  ) 

WRI  I  E  TAPE  2, ( Y(  I  )  ,  1=1 ,6  )  , X 

8  PRINT  14, RUN, ( P (  I  )  , 1=  1  , S  )  , ( YE (  I ) , YZ(  I ) , YR ( I)  ,  I  =  1  , 5 ) 

14  OFORYAK  3HRUN2X,  F4 .0//  1HPSX,  I  PSE  1  4 . 2/ /3HY/F3X »  IP. -SI  1  4 . 2/ /:>.  ;X/  3X, 

I  IP  3E  1 4 . 2//4i  IRLN  T2X ,  1 P  3F  1  4 . 2  /  /  3HS  *-*l  j3X  ,  1P3E  i  4 . 2  / /?HX  I  4  X  ,  1P3H4.2///) 
PRINT  80  ,  SPDfTX  ,  XF.SPU 

8  0  FOR.'-. AI  (6HMXSPD=2X,  1PE  1  4 . 2//o|  iXMSPl'J=2  X  ,  i  M-  14.2//) 

PRINT  40, ( Y(  I  )  t  1=  1 , 31 ,  ( Y (6 ) ) ,REV,  X  ,  Y ( 3  )  ,  Y ( 4 ) 

40  OF  ORE  AT  (  1  21  IF  I  N  AL  VALUFS//4HY/E  =  2X,  IPt  I  4. 2 //4HX/I>2X  ,  I  PE  14 . 2// 

1  6HTHET  A=2X  »  1  PE  1  4 . 2  /  /  4  H  S>J !  =2  X  ,  1  PL  1  4 .2  / /4I.KE  V  =  2X  ,  I  PE  14.2// 

2  2HX-2X,  IPE  I4.2//4HVEL  =  2X,  H’2F 14.9// ) 

EGYIN=2.*Y( 3) “P( 2)/P( 1 ) 

FGY  SUM= Y  (  6  )  *--2/P  (  1  )  +Y  (  2  )  **  2U  (  3  )  *Y(  1  )  **2+Y{ S  )  *«2  ►YU  )  **2 
1-2. »Y ( 6  )  *{ Y( b  )  »SINF ( Y(  => )  )-Y i 4 )*COSF(  \  ( 3 )  ) ) 

ER  ROR=  (  EGY  I  N-f  GY  SUE  )  «•  1  Ou  .  /  F  i.,Y  I  N 
PRINT  31  ,  EGY  IN,  F GY SUi  1 ,  ERROR 

3  1  OFORNAT  (6HEGYIN=  3X,  IFF  I  4 . 2  /  /  / 1  i  E  G  Y  >  U 2  X  ,  !  PF  14.2/  /  14HPE  RCENT  L'RK0R= 

1  2X,  IPE  14.2////) 

39  IF ( SENSE  SWI 1 CH  1  )  1  I  , 8 

3,.  I  F  (  F  ELP  )  1  1  ,  i  I  ,  12 

1  1  PAUSE  1 

READ  60, P (  I  )  , P ( 2  )  , DEL  P ,  Z  ,P  (  G  >  ,  P  (  4  ) 

60  FORi-  AT  (  6E8.2  ) 

GO  TO  SO 

12  IF ( Z-P ( 2 ) )  11 ,  1  1  , 30 
30  P(  2  )  =  l)ELP*P(  2  ) 

GO  TO  SO 


END 

SUBROUTINE  DEIUV(D,Y,X) 
DIMENSION  D(6)»Y(6),P(6) 


COMMON  P 
C  l  =  S  I N  F  I  Y  I  3  )  ) 
C2=C.OSF(  Y(  3)  ) 
C3=Y ( 6 ) «Y ( 6 ) 
C4=C2*  02 
U  (  1)  =  YI 4  ) 

D  (  2  )  -  Y  (  3  } 
D(3)=Y(6) 

0  D ( 4 ) - { - P (  I ) *C3 

1  «C2 

2  -P(6)*P(4 


3*C1*CI  +P( 

4C  \»C  1 
0  D I  3  )  =  I  C  2 
1C 1*C 1 

2*  P  I  3 ) »Y I  1 >*C1 
3  +  PI  I  )»P(4)*Y( 
4P(5)*C1 


+  P{  1 
)  »:  Y  (  4  )  -  P  ( 


C  1 

)*Y(2) 
)  *  Y  (  1 


0  0(6) 

1  *C2 

2  1  ./  (  I 


( P  t  2 ) »P I 


•PI 


-PI 

1  )  ) 


*04  -HMD* 

C  1  *C2  1-C..S 

-P(2)*Pt 5)*C?  -PI 

1  )  *P  (  6  )  »P  I  4  )  a  Y  I  4  )  *C  1  *  C  1  +  P  I  1  ) 

)*!./( l.-PI 1 ) ) 

«C 3  -P ( 4 ) *Y ( 5 )- Y I  2 ) -P (  1 ) *C2 

+  P  (  1  )  *  P  I  6  )  *  P  I  4  )  *  Y  (  4  )  *  C  I  -K-C2 

*C2  - P  (  1  )  *  C 2 «•  C 4 

) *C4  +PI  1  ) *  Y ( 2 ) *04 

)«•  l./l  l.-PI  1)1 
IfPI  1  )*P(6)*PI4)*Y(4)»C2 


*0  1 

)*C3 

PI  3 ) *Y  I  1 


-C  1 


)  » 


+  PI  1  ) 


*  c : 


+  PI  2) 


I  )*P ( 4 }*Y( 3) *C 1 


■PI  1  )  » Y 


+  P  I  1  ) 
I  2  )  »  C  1 


PI  3) 


■  Y  (  1  ) 
)  * 


END 

SUBROUT  IMF  RKUTTA I  NUMBER ,  XV  \R  ,  WARS, STEP ) 

DIMENSION  YVARS I  30) , AK( 4,30  )  ,  DYI 3 0 ) , YC I  30 )  , C I  3 ) 

C(  11=0.3 
C( 2  1=0.5 
C I  3  )  =  1 . 0 

CALL  DC  rVIDY,  YVAR  S » XV  A  P ) 

DO  1  J-  NUMBER 

1  AK ( 1 , J )  5TEP*DY(J ) 

DO  2  1=2,4 

XC=  XV  AR  +CI  1-  1  )*STtP 
DO  3  J=i, NUMBER 

3  YC  I  J  )  =  YVARSIJ)  ■»•  Cl  1-1  )*AK( 

CALL  DERIV(DYjYCjXC) 

DO  2  J=  1  ,  NU."  BE  R 

2  AKI  I  ,  J  )  =  STFP  «  D  Y I J ) 

DO  4  J=  1  , NUMBER 

4  YVARS  (  J  >  =  Y  V  A  R  S  (J)MAKI  1  ,  J  )  +2.*AK(  2,  J  )  +2  .  *  AK  I  3,  J  )  +AK  (  4  ,  J  )  1/6.0 
RETURN 

END 


• 

VALUES  or  P(6)  FOF.  P(3)  =  1.00 

PC  A)  -  0.00  (  UNOAMPpO) 

P(4)  -0.02  (OAMPEb) 

PC/) 

UNSUCC  ESS  PUL. 

SUCCESS  FU  L 

unsuccessful 

SUCCESSFUL. 

,  /.OOx/O'2 

0.00  x  IO~'2~ 

(e-3 0  x  /o'  z 

61 /  x /o'7- 

6.oo  x /o'2- 

3 .00  A  /O'3 

2.<k>Z 

2.14 

7A9 

2,0  2 

!  !‘O0 

bZG 

y.  3/ 

1,2-6 

1,31 

i  3.00  X/O '4 

6.1  S  A  10  ~ 5 

b.OO  x/O'3 

S.'X  X/O'3 

S-UU  x >o  ~3 

I,  00 

L - . - 

2.30 

2.9  2. 

2,2  T 

2,  4° 

I  3.  oo  X/O 

b2<i 

1.30 

8  CL  x  /o  ~4 

4 ,9o  x/o 

I  boo 

LoO  X /o'4 

0.30  x/0-4- 

4-0  3 

4,z4 

4> 


VALUES  of  pp 2)  FOR.  PC3)  ^0-26 

PCc)=  0-00  (UNOAMPEO) 

PC4 )-  0.02  (  DAMPEb  ) 

!  po ) 

UNSUCCESSFUL 

SUCCESSFUL 

UNSUCCESSFUL 

SUCCESSFUL 

1,00  x/o'7- 

3.1*8  x/o-'2- 

3-33  x  /O  "2 

3-62  X/O'2- 

3-  UP  x  / o'2- 

3.00  X  /O'3 

1,0  3 

hit 

/.EE 

1.03 

1,00 

<3,00  x  /o'3 

%-4o  x  /o'3 

1.2E  X  /0~3 

1,0/  x  /O'3 

3-00  X/O  ~4 

3-02. 

3.80 

3.  OS 

3-  20 

l,  00 

hi  3 

b  82 

J.  40 

h47 

3.00  x  /o  -p~ 

1-90  x/o  '4- 

?.2 o  x/o'4 

S.fco  x  JO  ~ ^ 

S.40  X  )0  '4 

i  i.oo 

3.13 

3.96 

2.  27 

2.  3% 

POINTS  DEFINING  TitE  BOUNDARY  t>LEWL' N  RECU..1  cF  SiwEGofoL 
AND  UNSUCCESSFUL  ACCELERATIONS 


\ 

values  of  pc?-)  fop  PCs)  -  o.  is 

PC4)  -  0,00  (  Cl  A/ DAM  pod) 

PC 4)  ^0.02.  (  DAMPSJU ) 

PO) 

UNSUCCESSFUL 

SUCCESSFUL 

unsuccessful 

SUCCESSFUL. 

/  .00  x  /O'2- 

4.17  x>D'z 

4-38  x  /o  ~z 

3.47  xio 

4-41  ?  10  "2 

?.oo  x  to '3 

4.04 

4.  SO 

3-lo 

4.04 

b-OO 

3.24 

3.4S 

3-/3 

3.28 

4- &o 

— 

— 

2-5f 

2-12 

4-30 

— 

— 

hi  4 

2.03 

4 .  oo 

0.38 

2. SO 

hll 

A  7  b 

3 -40 

2.08 

2.  tl 

— 

— 

3-00 

1-44 

bSb 

A  4!  X/o  "z 

l-4<4  x  io'7' 

mo 

— 

— 

?.?7  X  to  ~2 

1.33  x/o  "3 

/,  oo 

b.4Z  x/o  -3 

i ms  x  to ' 3 

S‘?7. 

b-/2 

3.00  a  /o'4- 

2-lb 

2%1 

2H0 

2.73 

A  OO 

A  27 

A  33 

A 40 

1-47 

3.  oo  jf  to 

s,io  x  /0'4- 

A-2 0  x  io  ~4 

sm  x/o  ~4 

b.oo  XJO'4 

2,to 

Co- 05 

4.30 

— 

- 

/■so 

4.1(0 

JTOO 

- 

— 

A  00 

3-7  8 

s,  IS 

2m  x/o  -4 

2.4  0  X/O'4- 

Points  Defining  the  Boundary  (Cent „ ) 


f 


VALUES  OF  PCI)  FOP  P(5 )  -  O.ro 

P(4)  -  0.  OO  (uNOA  MPFO ) 

/Y4)  =  A).  <02  £  DAMPED) 

pen 

U NS UCCESSFU  U 

SUCCESSFUL 

UNSUCCESSFUL 

SUCCESSFUL 

A  oo  X/O 2- 

A  9 S’  X  /0  ~z 

2.0s  x  /o  ~2 

1.76  X  /o  -2- 

b 75"  x  JO  ~2- 

7 -00  X /o~s 

/•54 

A  62 

2-6  7 

2.80 

S.Oo 

A  >9 

i,2<r 

A  <73 

7.08 

3.80 

9-1  /  x/0~3 

/■02 

1.41 

!.49 

3.00 

/■/I  x/o~z 

bn 

1.0/ 

A 

2.  (oO 

A  4-3 

mo 

A  27 

1.33 

2.20 

1,24 

1.30 

- 

— 

A  80 

/.OS 

l./O 

/.OS'  X/O  ~7- 

b/O  X  /O  -2' 

1-30 

9. 04.  x /o'3 

9.44  x  jo  '5 

- 

— 

1.20 

r.5-7 

9-oo 

- 

- 

1.00 

1-S4 

7. 90 

7-  /£>  ^  /O 

7.5^  X/O'* 

3.00  X  /0  ~4 

3*6/ 

3.80 

3.  !4 

330 

1,00 

A  72 

/.X/ 

A  42 

1.49 

5,00  x  to  ~ ^ 

7.6/  x/D'4 

8.00  x/0'4 

S-7/  X/0  '4 

6,-00  x /o'4 

I.  00 

3.72 

sno 

2-2  2 

2. 33 

2.&0  X  JO  ~3 

- 

- 

7.  <96  ^  /O  ~s 

9-S2  x  /o  ~3 

2.  30 

- 

- 

A  27  *  /o 

!■  33  x  Jo 2 

Points  Defining  tha  Boundary  (Cent*) 


APPENDIX  V 


MAXIMUM  AMPLITUDES  OF  VIBRATION  FOR  SUCCESSFUL 
ACCELERATIONS  THROUGH  THE  CRITICAL  SPEED  REGIONS 


1  P( 3)  -  J.  OO  P(4)  =  0-00 
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